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We propose a search for dark photons A' at the LHCb experiment using the charm meson decay 
-D*(2007)° —>■ D°A'. At nominal luminosity, decays will be produced at about 700kHz 

within the LHCb acceptance, yielding over 5 trillion such decays during Run 3 of the LHC. Replacing 
the photon with a kinetically-mixed dark photon, LHCb is then sensitive to dark photons that 
decay as A'—>■ e^e“. We pursue two search strategies in this paper. The displaced strategy takes 
advantage of the large Lorentz boost of the dark photon and the excellent vertex resolution of 
LHCb, yielding a nearly background-free search when the A' decay vertex is significantly displaced 
from the proton-proton primary vertex. The resonant strategy takes advantage of the large event 
rate for D*° —>■ D^A' and the excellent invariant mass resolution of LHCb, yielding a background- 
limited search that nevertheless covers a significant portion of the A! parameter space. Both search 
strategies rely on the planned upgrade to a triggerless-readout system at LHCb in Run 3, which will 
permit identification of low-momentum electron-positron pairs online during data taking. For dark 
photon masses below about 100 MeV, LHCb can explore nearly all of the dark photon parameter 
space between existing prompt-A' and beam-dump limits. 


I. INTRODUCTION 

Rare decays of mesons are a powerful probe of physics 
beyond the standard model (SM). Precise measurements 
of branching fractions and decay kinematics indirectly 
constrain extensions of the SM by bounding symmetry- 
violating or higher-dimensional operators. More directly, 
non-SM particles could be produced in meson decays 
when kinematically allowed, and depending on their life¬ 
times, these particles could yield striking signals with 
displaced vertices. A well-motivated hypothetical parti¬ 
cle is the dark photon A' which inherits a small coupling 
to the SM via kinetic mixing with the ordinary photon 
7 [IHSl- Indeed, some of the most stringent constraints on 
the properties of dark photons come from rare decays of 
mesons, including 7r°—7A' [7Ul4]. 7/77'—>■ 7A' [l5l[T6] . 
and <f) -T rjA' [T71 [TS] . 

The minimal dark photon scenario involves a single 
broken C/(l) gauge symmetry, along with mixing be¬ 
tween the A' and SM hypercharge fields via the oper¬ 
ator After electroweak symmetry breaking and 

diagonalizing the gauge boson kinetic terms, the dark 
photon gains a suppressed coupling to the ordinary elec¬ 
tromagnetic current where the relevant terms in the 
Lagrangian are 

C D + \m\,A'^A'^ + eeA/(1) 

This minimal scenario has two free parameters: the dark 
photon mass itia' and the kinetic-mixing parameter e (of¬ 
ten reported in terms of e^). The constraints placed on 
dark photons in the mA'-(^ plane are shown in Fig.[^ 
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for 2me < ruA' < 5 GeV, assuming that the A' domi¬ 
nantly decays into visible SM states (see Ref. (TH] for a 
review) 0 For > 10“®, the most stringent bounds come 
from searches for prompt A' decays at collider and fixed- 
target experiments [CT150115^ . As e decreases, the A' life¬ 
time increases, while as rriA' decreases, the lifetime and 
Lorentz boost factor both increase. Therefore, the con¬ 
straints obtained from beam-dump experiments exclude 
wedge-shaped regions in the rriA'-e^ plane [51 IT^ [TFl 1551 - 
[i5| . Also shown in Fig.[^ are electron g — 2 bounds 
[4511460 the preferred region to explain the muon g — 2 
anomaly m, and supernova bounds from cooling [55] 
and emissions |49j . Anticipated limits from other planned 
experiments are shown later in Fig.0 

In this paper, we propose a search for dark photons 
through the rare charm meson decay 

D*o^D°A', A'^e+e", (2) 

at the LHCb experiment during Run 3 of the LHC (sched¬ 
uled for 2021 -23)0 The goal of this search is to explore 
the region between the prompt-A' and beam-dump lim¬ 
its for the range G [10,100] MeV, which roughly in¬ 
cludes G [10“^°, 10“®]. Reaching such small values 
of is only possible for decays where the yield of the 
corresponding SM process (i.e. replacing A' with 7) is at 
least 0(10^®). Within the LHCb acceptance, over five 
trillion D*^ —> Z?®7 decays will be produced in proton- 
proton (pp) collisions at 14 TeV during Run 3, making 
this decay channel a suitable choice. 


^ There are also interesting searches where the dark photon decays 
invisibly to dark matter I20H29I . 

^ Since we follow the analysis in Ref. m, we obtain more conser¬ 
vative bounds from {g — 2)e than shown in Ref. IS]. 

^ Throughout this paper, = Z)*(2007)‘^ and the inclusion of 
charge-conjugate processes is implied. 
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FIG. 1. Current bounds on dark photons with visible decays 
to SM states, adapted and updated from Ref. m- The upper 
bounds are from prompt-j4' searches while the wedge-shaped 
bounds are from beam-dump searches and supernova consid¬ 
erations. The LHCb search region in Fig. [^covers most of the 
gap between these bounds for niA' < 100 MeV, with a reach 
extending to rriA' fS 140 MeV. Anticipated limits from other 
planned experiments are shown in Fig.[^ 

The range of niA' values that is in principle accessible 
in this search is mA' € [2TOe, Atod], where |50j 

ArriD = rrijjto — rrijjo = 142.12 ± 0.07 MeV. (3) 

The proximity of Amo to mT^o leads to phase-space sup¬ 
pression of the decay which results in a 

sizable branching fraction of about 38% for the decay 
D*o The small value of ArriD, however, also 

leads to typical electron momenta of ©(GeV) within the 
LHCb acceptance. Therefore, the planned upgrade to 
a triggerless-readout system employing real-time calibra¬ 
tion at LHCb in Run 3 m —which will permit identi¬ 
fication of relatively low-momentum e+e“ pairs online 
during data taking—will be crucial for carrying out this 
search. 

To cover the desired dark photon parameter space, we 
employ two different search strategies, shown in Fig.[^ 
The displaced search, relevant at smaller values of 
looks for an H'—)■ e+e” decay vertex that is significantly 
displaced from the pp collision. This search benefits from 
the sizable Lorentz boost factor of the produced dark 


^ This explains why we choose the decay D*^ —> D^A' instead of 
D*(2010)+^> D+A’, since the corresponding branching fraction 
^ D +7 is only 1.6%. 



FIG. 2. Potential bounds from LHGb after Run 3, for both 
the displaced (pre-module, solid blue) and resonant (dashed 
blue) searches. Also shown is an alternative displaced search 
strategy (post-module, dotted blue) that looks for A' vertices 
downstream of the first tracking module. 

photons and the excellent vertex resolution of LHCb. 
Our main displaced search looks for A' decays within 
the beam vacuum upstream of the first tracking mod¬ 
ule (i.e. pre-module), where the dominant background 
comes from misreconstructed prompt D*^ -A D^e~^e~ 
eventsj^ Because the A! gains a transverse momentum 
kick from pp collisions, the A! flight trajectory intersects 
the LHCb detector, making it possible to identify dis¬ 
placed e“'"e“ pairs with smaller opening angles than the 
HPS experiment [55] . We also present an alternative dis¬ 
placed search for A decays downstream of the first track¬ 
ing module (i.e. post-module), where the dominant back¬ 
ground comes from -A events with 7 -A e'^e~ 
conversion within the LHCb material. 

The resonant search, relevant at larger values of e^, 
looks for an A -A e+e” resonance peak over the con¬ 
tinuum SM background. This search benefits from the 
large yield of D*^ -A D^A decays during LHC Run 3, 
which is larger than the A yield in fixed-target exper¬ 
iments like MAMI/Al [SO] [3l] and APEX [32]. Fur¬ 
thermore, the narrow width of the D*^ meson, which 
is less than the detector invariant-mass resolution, pro¬ 
vides kinematical constraints that can be used to im¬ 
prove the resolution on m-g+g- . This resonant search can 
also be employed for non-minimal dark photon scenarios 


® We thank Natalia Toro for extensive discussions regarding this 
background. 
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where the A! might also decay invisibly into dark matter, 
shortening the AA lifetime. In that case, the anticipated 
limits in Fig.would roughly apply to the combination 
£2 X B{AA^ e^e-). 

The remainder of this paper is organized as follows. 
In Sec.[llj we estimate the A signal and SM background 
cross sections, extracting the production rate and 
decay modes from an event generator and estimat¬ 
ing the D*^ decay rates using a simple operator analysis. 
In Sec. m we describe the LHCb detector and charged- 
particle tracking, provide the selection requirements ap¬ 
plied to and D*^ meson candidates, and derive the 
A mass resolution. We present the pre-module displaced 
A search in Sec. IV a post-module variant in Sec.|^ and 
the resonant A search in Sec. VI Possible improvements 
are outlined in Sec. |VIl| and a comparison to other exper¬ 
iments (especially HPS) is given in Sec. VIII We sum¬ 
marize in Sec. |IX| and discuss how the LHCb dark pho¬ 
ton search strategy might be extended above the Amu 
threshold. 


II. SIGNAL AND BACKGROUND RATES 

Dark photon production in meson decays proceeds 
mainly via 


D*^^D^A, D%^{jA), (4) 

though the low-energy photon in the latter decay is un¬ 
likely to be detected at LHCb. Here and throughout, we 
use the notation X{YZ) to mean V —)■ YZ in a subse¬ 
quent decay. Because Amn < A —)• e+e“ is the 

only relevant visible decay channel. 

The dominant ba ckgr ounds to the pre-module dis¬ 
placed A search (Sec. IVI are i4*°—)■ D^e^e~ and D*^ ^ 
II° 7 r°( 7 e“'"e“), where the e“'"e“ pair is misreconstructed 
as being displaced due to a hard electron scatter in ma¬ 
terial. These backgrounds can be highly suppressed by 
requiring that the e+e“ kinematics are consistent with a 
displaced A vertex occurring in the proper decay plane. 
The dominant background to the post-module displaced 
A search (Sec.[^ is where the 7 converts 

into an e^e“ pair during interactions with the detector 
material. This background can be highly suppressed by 
requiring that the e“'"e“ vertex position is not consistent 
with the location of any detector material. The dominant 
backgrounds to the resonant search (Sec.|VI[) are again 
D*o £) 0 g-i-g- jj*o L)°7r°(7e+e“), where the 
A has been replaced by an off-shell 7 *. The first back¬ 
ground is irreducible, making the resolution on rue+e- 
the driving factor in the resonant search reach. 


A. D*^ Meson Production 

We simulate D*^ production in pp collisions at a 
center-of-mass energy of 14 TeV using Pythia 8.201 [S3] 


with the default settings. Since a large fraction of charm 
quarks are produced from gluon splitting and since we 
need to model forward physics at small transverse mo¬ 
mentum pt, we run all soft QCD processes in Pythia 
(i.e. SoftQCD;all = on). While the D*'^ production 
cross section is not yet known at 14 TeV, the result ob¬ 
tained using Pythia for the inclusive pp — l pro¬ 
duction cross section at 7 TeV agrees with the measured 
value by LHCb [51] to within about 5%j^ Since Pythia 
does not record the spin of the mesons, they are 
treated as unpolarized in this analysis. 

To define the fiducial region, we require the me¬ 
son to satisfy the following transverse momentum and 
pseudorapidity requirements: 

Pt{D°) > 1 GeV, 2 < r]{D°) < 5. (5) 


Note that this requirement is placed on the meson, 
not on the D*°, to suppress backgrounds to the com¬ 
ponent of the signal. The production cross section 
within this fiducial region is 

a{pp D*° Dgd) = 0.95 mb, ( 6 ) 


excluding secondary production of D*^ mesons from b- 
hadron decays. It may be possible to make use of some 
secondary decays; in this analysis, however, we require 
that the A originates from the pp collision to suppress 
backgrounds (see Sec. IV A). 

The nominal instantaneous luminosity expected at 
LHCb during Run 3 is 2 nb“^ per second [5T], which will 
produce D*'^ mesons at a rate of almost 2 MHz (equiva¬ 
lently, I II °7 at 0.7MHz). Assuming an integrated 
luminosity of 15 fb”^ in Run 30 this results in an esti¬ 
mated yield of 14 trillion mesons produced within 
this fiducial region, or 


7V(D*° ^ D° 7 ) = 5.4 X 10^2, (7) 


which we use as the baseline for our estimated reach. 


B. Meson Decays 

The meson is an I{J^) = |(1~) state with a mass 
of 2006.96 ±0.10 MeV and a width less than 2.1 MeV. It 
decays promptly mainly into two final states with branch¬ 
ing fractions of 

B{D*° D° 7 r°) = (61.9 ± 2.9)%, ( 8 ) 

B{D*^ D° 7 ) = (38.1 ± 2.9)%, (9) 


^ During the final preparation of this article, LHCb presented the 
first prompt charm cross section measurement at 13 TeV m- 
Based on this result, we estimate that the relevant cross sec¬ 
tion for determining the dark photon reach should be about 20% 
higher than the one used in this paper. 

^ The length of Run 3 is scheduled to be about the same as Run 1. 
LHCb collected a total of 3fb~^ in Run 1. The instantaneous 
luminosity will be five times higher in Run 3. Therefore, assum¬ 
ing the LHC performance is the same (including the slow ramp 
up), this gives an estimate of 15fb“^ in Run 3. 
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where the meson is a ^ (0 ) state [SD] . As mentioned 
above, _D*°—is the dominant background to the 
pre-module displaced search as well as to the resonant 
search. To our knowledge, this branching fraction has 
not yet been measured; therefore, we will estimate the 
rate for this decay using an operator analysis. This same 
approach is used to determine the —>■ Z3°A' rate. 

To calculate these D*^ —>■ transition amplitudes, 

we must first determine the (D*°| matrix ele¬ 

ment. By parity, time reversal, and Lorentz invariance, 
this transition dipole matrix element can be written in 
the form 


= Meff(fc') (10) 

where Va is the four velocity of the D*^ meson, kp is the 
momentum flowing out of the current, and e\ is the polar¬ 
ization of the meson. Here, /ies is a fc-dependent ef¬ 
fective dipole moment, whose value could be determined 
using a simple quark model (see, e.g.. Ref. [SS]) or us¬ 
ing a more sophisticated treatment with heavy meson 
chiral perturbation theory (see, e.g., Ref. m)- For our 
purposes, we simply need to treat /Teg as being roughly 
constant over the range G [0, Arnjj], which is a reason¬ 
able approximation given that Amjy < (Indeed, 

this relation is always satisfied in the heavy charm quark 
limit, where Amo c>c Aqqq/ me.) The precise value of 
/ieff is irrelevant for our analysis since it cancels out when 
taking ratios of p artial widths. 


Using Eq. ( |10[ ), we estimate the decay rate for D*^ — 
D^'y within the SM and in the Amu ^ mu limit to be 

OEM 




( 11 ) 


where oem = I^-k. To calculate the D*^ —>• 

decay rate, the off-shell photon propagator must be in¬ 
cluded. In the me = 0 limit, the amplitude for this pro¬ 
cess is 


\M 


2eVe'ff ' 


D*0-^D0e+e- \ = 


1 - 


{ki-vf+{k2-v) 


21 


ki ’ k2 


,( 12 ) 


where ki and k 2 are the electron and positron momenta. 
The ratio of partial widths is determined numerically to 
be 

D°e+e-) . , , 

= 6.4 X 10"^ (13) 

Since the dark photon also couples to 

Eq. (101 to calculate the L)*° —y D^A' decay rate. The 


ratio of partial widths is 
r(D*°^ D°A') 


r(£»*o^Ui07) ^ V Ami 


/ \3/2 




(14) 


where we assume m^', Am d ^ mu- This expression 
has the expected kinetic-mixing and phase-space suppres¬ 
sions. Since the D*^ meson is treated as unpolarized in 
Pythia, we ignore spin correlations in the subsequent 
A'—> e+e“ decayj^ 


As a technical note, to generate D*^ —^ D^A' events, we reweight 


C. Rare 7r° Decays 


To determine the —?► D^Tr^^-yA') decay rate in 

Eq. we start by estimating the rate of the decay 
7r° —)• jA' using the SM effective Lagrangian 

C = (15) 

Jtt 

where is the pion decay constant and the pion form 
factor is ignored. The dark photon is accounted for by 
making the replacement 

Ff^u —t Ff^i, + (16) 

which leads to the ratio of partial widths 


r(^°^7A0 

r( 7 r*^ -A 77 ) 


mZ 


ml 


3 


(17) 


The same effective Lagrangian can also be used for the 
SM decay 7r°—>• 76+e“. The amplitude is 


|A4 




4q?em 

Tr/^m^ _ 

•Z TT 'yp 


mlo + 2m^^_ -h m^+g- 


-b2m2 _m2+^_ -2mloiml^_ +ml+^-)]. (18) 


The ratio of partial widths is obtained numerically to be 


r(7r° -A 7e+e ) 
r(7r° -A 77 ) 


0.012, 


(19) 


which agrees with the nominal value for this ratio [50] . 


D. Dark Photon Decays 


Assuming the only allowed decay mode is A'—e“'"e , 
the total width of the A' is 


e^OEM 

La' = —— mA' 


m 


mt 


1 + 2-^ ! 1-4-^- (20) 


In the lab frame, the mean flight distance of the dark 
photon is approximately 

where 7boost is the Lorentz boost factor. In Fig.we 
show some example spectra of A' boost factors from sim¬ 
ulated —)■ D^A' decays, where both electrons are 


a sample of D*^ D®7 events from Pythia. In particular, we 

implement A' in the D*^ meson rest frame, boost to 

match the kinematics from Pythia, and then boost the D® 
decay products to account for the altered momentum. A 
similar strategy is employed for generating all other decays in 
our study. 
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14TeV: pp^ DL-4'(e+e-)„„ 



A. Track Types 

After exiting the VELO a distance of 0(1 m) from 
the pp collision, charged particles next traverse the first 
RICH detector (RICHl) before reaching a large-area 
silicon-strip detector located just upstream of a dipole 
magnet with a bending power of about 4 Tm |63j . Down¬ 
stream of the magnet, there are three stations of silicon- 
strip detectors and straw drift tubes. All tracking sys¬ 
tems will be upgraded for Run 3, though only the changes 
to the tracking systems upstream of the magnet are rele¬ 
vant here. The VELO has been redesigned to use pixels 
and is expected to have slightly better lifetime resolution 
and a lower material budget in Run 3 |64j . The tracking 
station just upstream of the magnet will also be replaced 
by a pixel detector and provide better coverage in p than 
the current detector [65] . This tracking station is known 
as the upstream tracker (UT). 

In LHCb jargon, there are two types of tracks relevant 
for this search: 


FIG. 3. Dark photon Lorentz boost factors for m^/ = 
{10, 20, 50,100} MeV. These factors are independent of e^. 


required to satisfy \pe\ = Pe > ^ GeV and 2 < pe < 5 so 


that they can be reconstructed by LHCb (see Sec. HI A 


below). The A inherits sizable momentum from the 
D*^ meson, leading to yboost factors that reach 0(10"*). 
The corresponding spectra for the total {£) and trans¬ 
verse (£t) flight distance of the A' are shown in Eig.|^ 
Eor < 10~^ the displacement between the A' and pp- 
collision vertices is resolvable by LHCb. 


III. BASELINE LHCb SELECTION 


• LONG tracks that have hits in the VELO, the UT, 
and the stations downstream of the magnet. These 
tracks have excellent momentum resolution in both 
magnitude and direction. 


• UP tracks that have hits in the VELO and the UT, 
but not in the stations downstream of the magnet. 
These tracks have excellent directional resolution 
obtained from the VELO. Since the curvature mea¬ 
surement is based only on the fringe held in which 
the UT operates, however, the uncertainty on the 
magnitude of the momentum is about 12% |65j . 


The LHCb detector is a single-arm spectrometer cov¬ 
ering the forward region of 2 < 77 < 5 jSS] [SS] ■ The detec¬ 
tor, which was built to study the decays of hadrons con¬ 
taining b and c quarks, includes a high-precision track¬ 
ing system capable of measuring charged-particle mo¬ 
menta with a resolution of about 0.5% in the region 
of interest for this search]^ The silicon-strip vertex lo¬ 
cator (VELO) that surrounds the pp interaction region 
measures heavy-havor hadron lifetimes with an uncer¬ 
tainty of about 50 fs m- Different types of particles are 
distinguished using information from two ring-imaging 
Cherenkov (RICH) detectors [6T|, a calorimeter system, 
and a system of muon chambers [62] . Both the momen¬ 
tum resolution and reconstruction efficiency are 0(10) 
times worse for neutral particles than for charged ones. 
For this reason, the analysis strategy outlined below is 
based entirely on charged-particle information. 


® The precision of electron momentum measurements is limited by 
bremsstrahlung radiation; see Sec.lIII D| 


We also note that LHCb defines DOWN tracks which 
have hits in the UT and downstream of the magnet but 
no hits in the VELO. While DOWN tracks are not used 
in this search, they could be useful for other searches 
involving long-lived particles. 

Charged particles may end up being reconstructed as 
UP tracks if they are swept out of the LHCb acceptance 
by the dipole magnetic field. This may occur if a particle 
is produced near the edge of the detector or if it is pro¬ 
duced with low momentum. For simplicity, we take any 
charged particle with 2 < p < 5 and p > 3 GeV to have 
100% efficiency of being reconstructed as a LONG track. 
Any track that is not LONG, but satisifies 2 < p < 5 
and p > 1 GeV is assigned as an UP track. In reality, the 
reconstruction efficiency is not a step function—particles 
with p < 3 GeV may be reconstructed as LONG tracks, 
while particles with p > 3 GeV may produce UP tracks 
or not be reconstructed at all—but this simple choice 
reproduces well the overall tracking performance. The 
momentum resolution for each track type is derived in 
App.0 







6 


14TeV: pp^ DgdA'(e+e-),„„ 14TeV: pp^ Dgj4l'(e+e-),„„ 




FIG. 4. Flight distance distributions for m^/ = {10,20,50,100} MeV showing (left) £ x (e^/10 ®) and (right) £t x (e^/10 ®). 


B. Reconstruction 

The D° meson momentum must be reconstructed for 
this search, since the kinematic constraints imposed by 
the D*^ mass will be used to suppress backgrounds and 
to improve the resolution on rrie+e-- We consider two 
categories of reconstruction. 

• F-type: All of the children are charged particles 
so that the can be fully reconstructed. At least 
two of the decay products must be reconstructed 
as LONG tracks. This suppresses combinatorial 
backgrounds and provides excellent resolution on 
the location of the decay vertex and on the 
momentum pjj. The remaining decay products are 
permitted to be reconstructed as either LONG or 
UP tracks. 

• P-type: At least two of the children are recon¬ 
structed as LONG tracks so that there is excellent 
resolution on the location of the decay vertex 
(there may be UP tracks as well). Requiring signifi¬ 
cant flight distance then permits reconstructing 
with good precision the direction of the momen¬ 
tum using the vector from the pp collision to the 
11° decay vertex. For the case where the invariant 
mass of the missing particle(s) is known, \pd\ can 
be solved for as discussed below. In this way, the 

is pseudo-fully reconstructed. 

The F^pe decays considered in this search are given 
in Table |l[ Each is of the form Zl° —>■ hh or 11° —>■ hhhh, 
where h = or tt^. We do not consider doubly 
Gabibbo-suppressed decays (e.g. 11° —>■ iF+7r“) since 


Decay 

B 

B X eff^ 

B X eff)) 

D°-> {K-tt+,KK, tttt} 

4.4% 

2.5% 

- 

{A-37r,2A27r,47r} 

9.1% 

4.5% 

1.0% 

Kl{v) 

6.8% 

- 

2.0% 

A^(^°)[o,™^ol 

22.0% 

- 

6.6% 

AA(A°)[,n] 

1.5% 

- 

0.5% 

A37r(^°)[o,™^„] 

8.5% 

- 

1.4% 

Total 


7.0% 

11.5% 


TABLE I. Decays of 11° mesons used in this search. The 
branching fraction B and efficiency-corrected branching frac¬ 
tion are given for each decay, for both the F-type (fully recon¬ 
structed) and P-type (pseudo-fully reconstructed) selections. 
The notation (a;) denotes that x is not reconstructed. Entries 
with an [o, 6] subscript count any decay where the invariant 
mass of the non-reconstructed system satisfies a < mmis < b 
as signal. 


they have small branching fractions and can be difficult 
experimentally to separate from the related Gabibbo- 
favored decays. LHCb has already published results 
using most of the F-type decays listed here (see, e.g., 
Refs. [Sni IHT]). and each decay is expected to have mini¬ 
mal combinatorial background contamination even with 
only loose selection criteria applied. Here we assume a 
baseline F-type Zl° selection efficiency of 90%. The total 
efficiency is then effp « 50%, which is dominated by the 
requirement that all decay products are reconstructed by 
LHGb. As shown in App.[Al the resolution on pjj for F- 
type decays is excellent. 
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In P-type decays, we can use the measured flight di¬ 
rection to pseudo-fully reconstruct po- The direction pu 
is a unit-normalized vector from the pp collision to the 
decay vertex. The magnitude is \pd\ = (Pvis’PD + 
Pmis'PD), where Pvis is the reconstructed (visible) mo¬ 
mentum and Pinis is the non-reconstructed (missing) mo¬ 
mentum. Balancing the momentum transverse to the di¬ 
rection of flight, requires = Ipmis x pd| = |pvis xpd\- 
Assuming that the invariant mass of the missing decay 
products is known, \pd\ can be solved for in the rest 
frame using conservation of energy and the known 
lueson mass. Since is invariant under boosts along 
Pd^ |Pmis’Pi>| in the rest frame is easily obtained. 
Finally, Pmis -PD can be determined in the lab frame up 
to a two-fold ambiguity that arises because the sign of 
Pmis’PD in the rest frame is not known. However, 
once the is combined with an A! —)■ e+e” candidate to 
form a D*'^ candidate, the vast majority of the time only 
the correct solution produces an invariant mass consis¬ 
tent with that of the D*^ meson. As described in App.[A| 
we take the baseline selection efficiency for P-type decays 
to be 50%, since the D'^ flight distance must be large rel¬ 
ative to the vertex resolution to obtain good resolution 
on pb- 

The P-type decays considered in this search are given 
in Table |Tj We again do not consider doubly Cabibbo- 
suppressed decays. Other decays that are ignored in¬ 
clude those where the missing mass cannot be reliably 
predicted, such as —>■ 7r~£~^{X), which dominantly 

has X = K^v as the missing system. Note that solv¬ 
ing for \pd\ in P-type decays requires using the known 
missing mass as a constraint. That said, the resolution 
is only degraded slightly if the true missing mass dif¬ 
fers from that used in the reconstruction by up to about 
0.2 to£)o- For example, when the visible part of the decay 
is , the most likely missing system is a single 7r°; if 

the missing mass is taken to be but the actual de¬ 
cay is K~'k^K^, the resolution obtained on We+e- 

by applying the “wrong” kinematical constraints to the 
D*° candidate is only worse by about 10%. In Table |lj 
we list the missing mass ranges considered as signal for 
each P-type decay. Candidates where the missing mass 
falls outside of these windows are ignored in this analysis, 
since they have worse resolution and anyways make up a 
small fraction of the P-type decays. A derivation of the 
P-type resolution is given in App.[A| The resolution 
on pb is about an order of magnitude worse in P-type 
than F-type decays; however, the mA' resolution after 
performing a mass-constrained fit is similar (as shown in 
Fig. below). 


C. D*^ Reconstruction 

To reduce the background from unassociated D^e~^e~ 
combinations, we require that the reconstructed mass dif¬ 
ference 

XrrAB° = TOreco(T’°e+e“) - mrecoiD°) (22) 


satisfies 


- 50 MeV < Amg™ - Attib < 20 MeV. (23) 

The looser requirement is placed on the lower edge due 
to bremsstrahlung by the electrons. This mass require¬ 
ment highly suppresses the decay D*^ ^ D^TT^{'ye~^e~) 
and its A! counterpart, except when me+e- is large (see 
Sec. IV D] below). The efficiency of this requirement is 
about effAmn ~ 85%. Note that this cut can be tight¬ 
ened at the expense of signal efficiency if combinatorial 


backgrounds turn out to be problematic (see Sec. IV C 
below). 


D. A! Reconstruction 


The reconstructed electrons produced in A! —>• e+e“ 
decays are a mixture of UP and LONG tracks. Only a 
few percent of the electrons have momenta large enough 
that equivalent-momenta non-electrons would be able to 
emit Cherenkov light in RICHl. Therefore, identification 
of the and e~ should be highly efficient with a low 
hadron-misidentification rate. Furthermore, the signa¬ 
ture of a maximum-Cherenkov-angle ring in coincidence 
with a track should suppress the fake-track background 
which can be sizable at low momenta. 

Bremsstrahlung radiation and multiple scattering of 
the electrons significantly affect the rue+e- resolution. 
We implement this numerically in our simulation fol¬ 
lowing Refs. [50l |68] and using the Run 3 LHCb 
VELO [M], RICHl [55], and UT |5S| material budgets. 
Bremsstrahlung downstream of the magnet does not af¬ 
fect the momentum measurement and is ignored. 

In Fig.[^ we show the resolution on rue+e- for several 
values of tua'-, where the A! candidates are constrained 
to originate from the pp collision. Bremsstrahlung cre¬ 
ates large low-mass tails resulting in poor resolution on 
me+e-- Since the mass is known and its width is 
less than the detector resolution, though, we can correct 
the me+e- distribution once we identify the candi¬ 

date and apply the ArrAff^° cut. As a heuristic, one can 
rescale the mg+e- value by a simple correction factor 

/ Am*'®™ \ 

mX-=<Te°-(2-^). (24) 


A more sophisticated approach involves performing a 
mass-constrained fit to enforce energy-momentum con¬ 
servation and the known mass using the covariance 
matrices of all reconstructed particles. Using this fit, we 
find 10-20% improvement in (j{m^+^-) relative to the 
simple correction given in Eq. (24). As shown in Fig.j^ 
the resolution on mg+e- after the applying the kinematic 
fit is 2-3 MeV using F-type candidates, and 2-5 MeV 
using P-type candidates. 

The key difference between the pre-module displaced, 
post-module displaced, and resonant searches are the re¬ 
quirements placed on the A' flight distance. These are 
described in more detail in the subsequent sections. 








14TeV: pp^ DL-4'(e+e-)„„ 



FIG. 5. Distribution of Tn^+^- with (solid, dashed) and 
without (dotted) incorporating the D*^ mass constraint for 
’’nA' = {10, 20, 50,100} MeV. The solid curve shows better 
performance than the dashed one because F-type candi¬ 
dates have better momentum resolution than P-type ones. 

IV. DISPLACED A! SEARCH (PRE-MODULE) 

The A! typically has a large Lorentz boost factor, re¬ 
sulting in the A! decay vertex being significantly dis¬ 
placed from the pp collision for < 10“^. The com¬ 
bined signature of a displaced decay vertex, a dis¬ 
placed A! —e+e“ vertex, m{D^A') consistent with 
m{D*^), and a consistent decay topology will result in 
a nearly background-free search. This pre-module dis¬ 
placed search is aimed at A' decay vertices that occur 
within the beam vacuum upstream of the hrst VETO 
module intersected by the A! trajectory. 

A. Conversion and Misreconstruction Backgrounds 

At LHCb, the first layer of material is the foil that sep¬ 
arates the beam vacuum from the VETO vacuum. This 
foil is corrugated to accommodate the VETO modules, 
such that if the AA decays prior to the foil, it still ef¬ 
fectively decays within the VETO tracking volume. The 
average transverse distance that the AA will travel before 
hitting a VETO module is 6 mm [53], which, because 
of the corrugated foil geometry, is roughly the average 
transverse flight distance to the foil as well. 

To effectively eliminate backgrounds from 7 —>■ e“'"e“ 
conversions in the foil, we require the AA decay vertex 
to be reconstructed upstream of the foil. Eurthermore, 
each reconstructed electron must have an associated hit 
in the first relevant VETO module given the location of 


the reconstructed AA decay vertex. These hits are re¬ 
quired to have at least one vacant VETO pixel between 
them to avoid any charge-sharing issues, imposing an ef¬ 
fective buffer distance between the A decay vertex and 
the foil: 

0.123 mm , , 

-, (25) 

e~ 

where cte+e- is the electron-positron opening angle. In 
reality, the VETO pixels in Run 3 will be 55 x 55 pm^ 
squares; the definition of D is based on treating the pix¬ 
els as circles with 0.123 mm being twice the effective di¬ 
ameter (the precise value used here has no impact on our 
search). The pre-module A' requirement can then be ap¬ 
proximated by requiring the A' transverse flight distance 
to satisfy 

It < 6 mm — Dx, L^x = D sin 9, (26) 

where 9 gives the A flight direction. To remove A tra¬ 
jectories that first intersect the foil far from a module, 
we require t\a' > 2.6. We also impose r\A' < 5 to avoid 
possible contamination due to pp collisions that are not 
properly reconstructed [^ 

Having suppressed conversion backgrounds, the dom¬ 
inant background comes from prompt D*^ —> D^e~^e~ 
events where the e“'"e“ vertex is misreconstructed as be¬ 
ing displaced because of multiple scattering of the elec¬ 
trons in the detector material. We estimate this back¬ 
ground in a toy simulation of the Run 3 VETO, taking 
scattering angle distributions from a Geant simulation 
which includes non-Gaussian Moliere scattering tails PI 
Many of these fake A vertices can be eliminated by re¬ 
quiring a consistent decay topology, in particular that 
the angle between pA' and the vector formed from the pp 
collision to the A decay vertex is consistent with zero, 
and the electrons travel within a consistent decay plane. 

The remaining misreconstructed background events 
have a consistent topology, so a cut on transverse flight 
distance £x is required to ensure a significant displaced 
A vertex. To avoid fake displaced vertices from one 
electron experiencing a large-angle scattering, we also re¬ 
quire both the electron and positron to have a non-trivial 
impact parameter (IP) with respect to the pp collision. 
These requirements are summarized by 

Ti 

£t > IPe± > (27) 


An A' candidate may be accidentally formed from a prompt 
e'*"e~ pair produced in a pp collision if the event is not prop¬ 
erly reconstructed. In particular, if a meson is produced in 
another pp collision upstream of that interaction point, the “dis¬ 
placed” A' would produce a consistent decay topology, albeit 
with Tjj^/ 00 . 

It is likely that Geant overestimates the probability for large- 
angle scatterings (see Ref. ESI). If SO, our results are conserva¬ 
tive, since these scattering tails effectively define the reach for 
the pre-module A' search. 
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FIG. 6. Pre-module displaced search A' signal efficiency 
effA'(^'A', £^) for several values. The dashed line shows 
the efficiency of the track-reconstruction requirements placed 
on the electrons alone, while the solid lines show the total A' 
efficiency. This efficiency does not include the contribution 
from the or selection (i.e. eff^’^ or effAmo). 


where the value of n is adjusted to yield « 1 background 
event in each A' mass window, with n ranging from 3 
to 5 as a function of toa'. The selection in Eq. (27) is 
meant to be simple and robust, and could certainly be 
optimized in a full analysis. See App.|^for details on the 
and IP resolution. 


B. Event Selection 

Summarizing, the event selection for the pre-module 
displaced A' search is: 

• F-type or P-type candidate; 

• and e~ from LONG or UP tracks with hits in 
the first VELO module they intersect; 

• reconstructed D*^ candidate from the 0°, e+, and 

e"; 

• reconstructed A' —>■ e^e“ satisfying the conversion 
veto < 6 mm — Dt, ^a' G [2.6, 5]); 

• reconstructed A' -A e+e“ with significant displace¬ 
ment (^T > n(Tf.j,,IPe± > fcip). 

In Fig.|^ we show the resulting A' signal efficiency 
effA'(mA',e^). For = 10“^°, which is near the low end 
of the reach, the efficiency is limited by the efficiency of 
the conversion veto. As e increases, the requirement of 


a significant displacement ultimately limits the reach of 
the displaced search. 


C. Additional Backgrounds 

Beyond the misreconstructed D*^ -A D^e^e~ back¬ 
ground, a full accounting of the potential backgrounds 
for the displaced search is difficult since all SM processes 
with large rates are highly suppressed. Therefore, any 
additional backgrounds will be dominated by extremely 
rare processes or highly unlikely coincidences. 

One possible source of backgrounds would he B ^ 
decays, since the resulting e+e“ vertex 
is truly displaced. Such decays can be suppressed by 
making the following requirements: the and A' mo¬ 
menta must intersect the pp-collision point when traced 
upstream from their respective decay vertices; the D*^ 
decay vertex is consistent with the pp-collision vertex; 
and there are no additional tracks consistent with origi¬ 
nating from the A' decay vertex. Furthermore, one could 
require that the A' decay vertex is downstream of the 
decay vertex, which would be efficient for the smaller 
e values probed in this search. Therefore, we do not 
expect a significant amount of background coming from 
B—)• ZI*°(L>°e+e“)A' decays. 

The decays of other long-lived mesons could also be 
sources of displaced e+e“ vertices. Decays of charged 
pions and kaons that produce an e+e” pair are rare, 
though, and the probability for these particles to decay 
in the VELO is small. A more likely source is the decay 
7r° -A- e“''e“7, where the 7r° is produced in the decay of 
a long-lived meson. All of the other meson-decay prod¬ 
ucts must be neutral, of course, otherwise the presence 
of additional charged particles consistent with originat¬ 
ing from the e+e“ vertex could be used as a veto. For 
example, the decays ATg—>■ 7r°7r° and —>■ AgTr® occur 
with huge rates within the LHCb VELO. Such decays, 
however, are unlikely to result in the A' candidate mo¬ 
mentum intersecting the pp-collision point or to occur 
in coincidence with a meson such that m{D^A') is 
consistent with m{D*^). 

To see whether we could estimate displaced combinato¬ 
rial backgrounds in Monte Carlo, we generated a sample 
of 30 million Pythia pp collisions at 14 TeV. We found 
that no combination of a true with two displaced 
tracks (not necessarily electrons, but assigned the elec¬ 
tron mass) had an invariant mass within the D*^ mass 
window. In this simulated sample, there are only three 
candidates with m{D^e'^e~) — m{D^) < 500 MeV us¬ 
ing true electrons and none within 150 MeV of our D*^ 
mass window. However, we are anticipating 10® times 
more meson decays in the full LHCb data sample, 
so it appears that it is not feasible to use Monte Carlo 
to precisely estimate the displaced combinatorial back¬ 
ground. This type of background will therefore need to 
be examined in data using the sidebands. If spe¬ 

cific sources of combinatorial background are identified as 
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problematic, then the selection will need to be adjusted 
to remove them. 


D. Contribution from Pion Decays 


Thus far, we have ignored the channel D*^ —>■ 
D^Tr^{'^A'), which is another potential source of sig¬ 
nal events. Decays of this type are highly suppressed 
by the requirement in Eq. (23), though, unless 


niA' is large. For most allowed niA' values, one can 
choose whether or not to include such decays in the 
analysis by adjusting the requirement. After re¬ 

moving the Arrf£'^° requirement, the expected yields of 
D*o £)*o £)0j^i (j0(jays are com¬ 

parable, but so is the expected background contamina¬ 
tion from misreconstructed D*^ —>■ Zl°7r°(7e“'"e“) and 
D*o jjOg+g- choose not to include this chan¬ 

nel when estimating the reach below, but note that such 
decays may prove useful in a complete analysis. 

If one does try to use the —>■ chan¬ 

nel, then one should be aware of an important subtlety 
when incorporating We+e- information. As described in 


Sec. Ill D a kinematic fit can be used to improve the 
rue+e- resolution. For D*^ —>• II°7r°(A'7) decays, how¬ 
ever, the missing 7 is not accounted for when enforcing 
energy-momentum conservation. We find that this re¬ 
sults in the vriA' peak being shifted up in mass by about 
20 MeV, with the resolution on mg+e- degraded by about 
a factor of two. This results in two peaks in the recon¬ 
structed mA' spectrum, coming from the D*'^ —> D^A' 
and )■ II°7r°(7A') channels. If the background level 
is low, then this second peak could be used to boost the 
significance of an A' signal. Indeed, one narrow peak with 
a second, wider peak shifted in mass by a fixed amount 
would be a striking signature. If the background level 
is high, then this wider peak would largely be absorbed 
into the background and have no impact on the signal 
significance. 

Finally, photon conversions arising from —> 

D°7r°(77) decays are also highly suppressed by the 
Arrf£^° requirement, and in the absence of misreconstruc- 
tion, can be eliminated by the pre-module requirement in 
Eq. ( [M| ). We expect such conversions to contribute less 
than those from —>■ T>°7 decays. 


E. Reach 


The expected signal yield for the displaced search as a 
function of ttia' and is given by 




D°A') 


~ 85 


10-10 


X (^eff^ +effo) effA'(mA',e^) 

3/2 


effAr; 


1 - 


Amjj 


eSA'(mA',t)- (28) 


As discussed above, we adjusted the requirement in 
Eq. (27) to ensure « 1 background event in any given 
niA' window. Assuming that all relevant backgrounds 
have been accounted for, the reach would be set at 95% 
confidence level by requiring S > 3. However, to allow 
for the possibility of extremely rare background sources, 
the reach is set by requiring S > 5. In this way, we ac¬ 
count for either additional background candidates in the 
final data sample or for a lower A! selection efficiency 
due to the criteria required to suppress these additional 
backgrounds. The reach is shown in Fig.[^ assuming 
15fb“^ of data collected by LHCb in Run 3, which cov¬ 
ers a significant part of the allowed parameter space for 
mA' ^100 MeV. 


V. DISPLACED A' SEARCH (POST-MODULE) 


In order to capture more A' signal events, one can ef¬ 
fectively reverse the pre-module requirements in Eq. (26) 


and search for post-module A! decays. Here, the dom¬ 
inant background is —>■ where the on-shell 7 

converts into e“'"e“ via interactions with the detector ma¬ 
terial. As we will see, this post-module search does not 
cover much additional A! parameter space compared to 
the pre-module search, but is important as a cross check 
of a possible AA discovery. 


A. Misreconstruction and Conversion Backgrounds 

The background considerations in the post-module 
case are reversed compared to the pre-module case in 
Sec. jlV A} Here, the background from misreconstructed 
D*o — D^e^e~ events can be effectively eliminated by 
requiring no hits in the first VETO module intersected 
by the reconstructed electron trajectories. 

The dominant background in the post-module search 
comes from D*^ -A D^y with photon conversions. We 
simulate this background using the Run 3 LHCb VETO 
material as described in Ref. (64] with the Bethe-Heitler 
rUe+e- spectrum as given in Ref. mE We start with 
electron tracks that each have at least three hits in the 
VETO. This imposes an effective A! transverse flight dis¬ 
tance requirement of 

G [6mm,22mm]. (29) 

In reality, the electron hit requirement does not result 
in a step function for the AA efficiency. However, in the 
long-lifetime limit, this simple approximation produces 


It is vital that all A! searches use the Bethe-Heitler spectrum, 
rather than the one produced by Geant. Geant vastly under¬ 
estimates the fraction of conversions that produce large 
due to the usage of a less-CPU-intensive approximation of the 
Bethe-Heitler equation. 
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FIG. 7. Post-module displaced search A! signal efficiency 
effA'(^'A', £^) for several values. The dashed line shows 
the efficiency of the track-reconstruction requirements placed 
on the electrons alone, while the solid lines show the total A! 
efficiency. This efficiency does not include the contribution 
from the or selection (i.e. eff^’^ or effAmo). 


the same integrated efficiency. We then require the re¬ 
constructed A! vertex to be significantly displaced from 
the VELO material. This can be well-approximated by 
treating the VELO as a stack of modules located at lon¬ 
gitudinal distances Zi = f • 25 mm, where z is measured 
from the point where the d' has = 6 mm (i.e. the 
average position where the A trajectory crosses a VELO 
module). From a given A decay vertex at a location z 
between modules i and i + 1, one requires 


z — Zi > naz, 




— z > D, 


IPi,e± > 


where D is the same buffer distance in Eq. ( |25[ ) and IP^ 
is defined with respect to the location where the A tra¬ 
jectory intersects the f-th module (see App.jB] for the 
corresponding resolutions). We also impose the same 
T^A' G [2.6, 5] requirement as in Sec. IV A Using our sim¬ 
ulation, we adjust n such that « 1 > ZI°7,7—)■ e+e“ 

event will survive these criteria in each A mass window. 


B. Event Selection and Reach 

Summarizing, the event selection for the post-module 
displaced A' search is: 

• F-type or P-type candidate; 

• and e~ from LONG or UP tracks with no hits 
in the first VELO module they intersect; 


reconstructed candidate from the e'*', and 


• reconstructed A —>■ e+e satisfying the prompt 
veto {£t G [6 mm, 22 mm], r]A' G [2.6,5]); 

• reconstructed A —>■ e+e” with significant displace¬ 
ment from VELO modules (Azi > ncz, IPi^e± > 
fo-jp., Azi+i > D). 

In Fig.[^ we show the resulting A signal efficiency 
effA' (wAqe^). 

The formula for the reach for the post-module dis¬ 


placed search is identical to Eq. (28), albeit with a modi¬ 


fied signal efficiency. We again set the reach using S' > 5 
to allow for rare, unaccounted for background sources. 
As shown in Fig.[^ the reach in the post-module search 
is not any better than the pre-module search. The rea¬ 
son is that for small enough e^, the decay probability 
as a function of is constant, and there is comparable 
efficiency for the A decay vertex to be upstream or down¬ 
stream of the first VELO module. As increases, the 
post-module £t > 6 mm requirement becomes inefficient 
before the pre-module vertex requirements do, giving the 
pre-module search a better reach. Of course, a slightly 
better limit could be set by combining the pre-module 
and post-module requirements, but we do not do this 
here since the dominant background sources are differ¬ 
ent. We also note that if a discovery is made, the ability 
to confirm the presence of a signal in both displaced re¬ 
gions will provide a powerful systematic check. 


C. Additional Backgrounds 

The same additional backgrounds from Sec. |IV Cj might 
affect the post-module search, with the exception of 7r° 
mesons produced in charm-meson decays. In addition, 
there is a potential background from improper recon¬ 
struction of photon conversion events. 

One way to misreconstruct a 7 conversion as an A de¬ 
cay is if the 7 converts in or just in front of a VELO mod¬ 
ule. Since the separation between the and e“ when 
they traverse the VELO module would be less than the 
hit resolution, only a single hit would be recorded. The 
positron track may be formed using this hit. The electron 
track would be missing a hit in this module and fail the 
selection requirements, unless an unassociated hit hap¬ 
pens to occur close by. The hit occupancy expected in 
the VELO during Run 3 is about 0.08% in the inner¬ 
most pixels, and less than 0.01% for pixels more than 
10 mm from the beam line [64) . Since we require a full 
vacant pixel between the first e"*" and e~ hits, there must 
be a large angle between the electron momentum and 
the vector h formed from the unassociated hit to the e~ 
hit in the second module. Therefore, a consistent track 
will only be formed if the electron undergoes an unlikely 
scatter in the second module so that the hit in the third 
module is consistent with h. We cannot reliably estimate 
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the probability of this coincidence in our toy simulation, 
but expect it to be at worst comparable to the remain¬ 
ing conversion backgrounds left after the requirements in 
Sec.|V A[ so it should not affect the predicted reach. 


VI. RESONANT SEARCH 


When is large enough, the A' decay vertex is no 
longer significantly displaced from the pp collision. In 
this case, we have to rely on reconstructing an A! mass 
peak. This search is also relevant for non-minimal sce¬ 
narios with a larger dark photon width, such as when 
the dark photon has an invisible decay to dark matter 
particles (see, e.g.. Refs. POMS])- 

For simplicity, we perform our analysis still assuming 
the minimal Fvalue in Eq. (201. This means that con¬ 
version backgrounds are important, so we apply the con¬ 
version veto from Sec. IV A The dominant background 
becomes prompt —>■ which is irreducible, 

and the resonant reach is limited by the rUe+e- resolution 
shown in Fig.[^ To ensure good rUe+e- resolution, we 
further impose an A' opening angle cut tte+e- > Smrad. 

In summary, the event selection criteria for the reso¬ 
nant search are 


• F-type or P-type candidate; 

• and e~ from LONG or UP tracks with hits in 
the first VETO module they intersect; 

• reconstructed D*^ candidate from the 0°, e+, and 

e"; 


• reconstructed A' —> e“'"e satisfying conversion 
veto (t^T < 6 mm — Ut, VA' G [2.6, 5]); 


• A decay opening angle > 3mrad to ensure good 
rUe+e- resolution. 

In Fig.[^ we show the A' signal efficiency eff^'(TO^',e^) 
for the resonant search. Apart from the electron track- 
reconstruction, the dominant source of inefficiency is the 
opening angle requirement. 

Due to the large irreducible background level in this 
search, we assume that —>■ D^'k^{'^A') decays pro¬ 

vide negligible additional sensitivity and therefore ig¬ 
nore them. Background contamination from D*^ -A 
II°7r°(7e+e“) decays is still included, but it is highly sup¬ 
pressed by the AnA£‘^° requirement except when m^+g- 
is large. 

To determine the reach in the resonant search, we re¬ 
quire S/y/B > 2 to obtain a 95% confidence limit. The 
signal significance is given by 


S 

71 


= D'^A') 


I N {D*o-)■ DO-f) 


' (effl,effA'(mA',e2), (31) 


D*^ ^ A'(e'^e ) [resonant] 



FIG. 8. Resonant search A signal efficiency eff^'(m^', e^). 
The dashed line shows the efficiency to reconstrnct both elec¬ 
trons, while the solid line shows the total A efficiency. The 
efficiency is nearly independent of e for > 10“^. This effi¬ 
ciency does not include the contribution from the or 
selection (i.e. eff^’^ or elf Amo). 


where 


AF = 


rmj^i +Am^i 


dr{D*°^ D°e+e-) 

dirip 


(32) 


and Am A' = 2 a{mg+e-) is evaluated at mA' 0 Here, 
we are assuming that eSA’imA'^A) is constant over 
[mA' — AmA>,mA' + AmA']- The reach for the resonant 
search is shown in Fig.[^ and combined with the pre¬ 
module displaced search, it closes most of the available 
A parameter space for mA' ^ 100 MeV. 

The triggerless-readout system, along with real-time 
data calibration, will make it possible to identify D*'^ -A 
D^A candidates online during data taking. For the dis¬ 
placed searches, the number of candidates recorded for 
further analysis can be made as small as required since 
the dominant backgrounds are reducible. The resonant 
search, on the other hand, relies on looking for an m^+e- 
peak in 0(10®) D*^ -A D^e'^e~ decays. Recording bil¬ 
lions of full events for such a search is simply not feasi¬ 
ble. That said, LHCb is already commissioning in Run 2 
partial event storage for use in studying high-rate pro¬ 
cesses [55]. This involves storing only the information 


Because F-type and P-type D® mesons yield different A! 
mass resolutions, Eq. | |32| l is evaluated using the appropriately 
weighted average of the two event categories. 
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relevant to the signal candidate and discarding the rest 
of the event information. 

In terms of required bandwidth, in principle the reso¬ 
nant search can be carried out by storing just rrie+e- for 
each signal candidate, though in practice one wants to 
store additional information to perform cross checks. At 
minimum, we expect that the four-momenta and various 
detector-response information can be kept for all parti¬ 
cles in each signal candidate. More information can be 
kept for displaced candidates that are inconsistent with 
arising due to photon conversion, where the criteria used 
to select such candidates will be much looser than those 
applied to define our reach above. Finally, full events can 
be kept for a small fraction of signal candidates to per¬ 
mit detailed offline studies of the detector performance 
and background contributions. The experience gained by 
LHCb during Run 2 using reduced event storage will aid 
in determining how to optimize the data-storage strategy 
employed for the D*^ D^A' search. 


VII. POSSIBLE IMPROVEMENTS 


The predicted reach in Fig. [^covers most of the avail¬ 
able A' parameter space for ttia' ^ 100 MeV, and one 
expects that LHCb could fully cover this region with 
modest improvements. That said, this reach is based on 
simulation, so despite our conservative approach, the ac¬ 
tual reach may be less than predicted. Furthermore, the 
ultimate goal is to discover an A! boson and measure its 
properties. Discovery requires 5 (t sensitivity, which our 
search provides over most of the relevant A! parameter 
space but not all of it. Precision measurement of the A! 
mass and lifetime may be vital to determining whether 
the A! decays invisibly into dark matter. For these rea¬ 
sons, it is worth investigating possible improvements. 

There are various ways to improve the dark photon 
search at LHCb. At minimum, one could assume that 
the integrated luminosity accumulated by LHCb is larger 
than the 15 fb“^ baseline used in this studyj^ One year 
of running at the nominal instantaneous luminosity cor¬ 
responds to about 10fb“^; therefore, our Run 3 esti¬ 
mate conservatively provides substantial LHC ramp up 
and LHCb commissioning time. LHCb also plans to take 
data in Run 4, collecting a minimum of 50fb“^ in Runs 
3 and 4 combined m- Furthermore, LHCb is currently 
investigating ways to increase the integrated luminosity 
in Run 4. By the end of Run 4, LHCb may collect 10-20 
times more data than we used to estimate the reach in 
this study. 

There are several ways to improve the We+e- resolu¬ 


tion. The mass-constrained fit we employed in Sec. HI D 


only makes use of kinematical information. For large 
flight distance, the vector formed from the pp collision 
to the A decay vertex provides another measurement 
of the direction of flight. Thus, a more sophisticated 
fitting procedure would not only constrain energy and 
momentum, but also enforce a consistent decay topol¬ 
ogy IZ2], likely improving the resolution. In addition, 
during Run 1 LHCb managed to recover many of the 
bremsstrahlung photons emitted by electrons by recon¬ 
structing them in the calorimeter system [73]. These 
photons were then added back to the electron momenta. 
Since the calorimeter occupancy will be higher in Run 3, 
we chose to assume that this recovery procedure would 
not be possible in our analysis. If bremsstrahlung recov¬ 
ery is possible, though, it would improve the rUe+e- reso¬ 
lution by 0(10%). We note that if a major upgrade of the 
calorimeter system |74| is installed for Run 4 that utilizes 
large-area picosecond photodetectors m. this would sig¬ 
nificantly improve the mg+g- resolution. Such an up¬ 
grade would also permit using many more decays. 

Installing additional tracking stations onto the face of 
the dipole magnet would provide LONG-track-level mo¬ 
mentum resolution for UP tracks m- The hit resolution 
required is ©(mm), which means that such a tracking 
station would not be too costly. Since many A decays 
are reconstructed with at least one UP track, installing 
such a system could greatly improve the rUg+g- resolu¬ 
tion. Furthermore, many decays produce UP tracks 
whose resolution would also be improved by installing a 
tracking system onto the face of the magnet. 

For simplicity, our resonant search imposed a conver¬ 
sion veto, but since the conversion background is smooth, 
one could still perform a resonant search without any 
requirement at the expense of larger backgrounds. This 
is particularly relevant for small mA' where the buffer 
requirement is most inefficient. 

Finally, one could consistently combine both displaced 
searches and the resonant search following a strategy sim¬ 
ilar to Ref. m- This is particularly relevant if the A 
can also decay into dark matter, since then t he r elation- 
ship between F^q toas and given in Eq. (20) will no 
longer hold and A will have a shorter lifetime. To avoid 
introducing model dependence into the A search, one 
can define exclusive prompt and displaced regions exper¬ 
imentally, and then optimally combine the information 
from both regions regardless of the A lifetime m- A 
combined strategy might also help close the remaining 
gap between the resonant and displaced searches at large 
niA'- 


VIII. COMPARISON TO OTHER 
EXPERIMENTS 


As stated in footnote [^ the first measurement of the cc cross 
section at 13 TeV suggests that the cross section in Eq. is 
likely about 20% too small. If so, one can view our reach results 
as for 12fb~^. 


There is a rich planned program of dark boson searches 
m, and various experimental proposals are sensitive to 
the same dark photon mass and coupling range targeted 
by this —>■ D'^A search. Their anticipated reaches are 
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FIG. 9. Comparing the LHCb reach to other proposed dark 
photon experiments. 

summarized in Fig.[^ Experiments like APEX [321 [77], 
MESA/MAMI [75], DarkLight [7ni]5n], VEPP-3 [13], and 
Mu3e |81j are high luminosity experiments that use a 
resonant search strategy]^ Experiments like SHiP [52] 
and SeaQuest |83j are beam-dump experiments that use a 
displaced strategy. See App.[^for discussion of a possible 
DO A' search at Belle-II. 

The closest comparison to the LHCb D*o —j, D^A' 
search is the HPS experiment [52]. Because HPS has 
a dedicated tracking and vertexing detector, it is able to 
search for both resonant and displaced A' signals, corre¬ 
sponding to the upper and lower HPS regions in Fig.]^ In 
terms of vertex performance and mass resolution, HPS is 
quite comparable to LHCbj^ It is therefore reasonable to 
ask why the resonant and displaced search regions over¬ 
lap for LHCb but not for HPS. 

There are three main advantages of LHCb over HPS, 
at least for tua' < 100 MeV. 

• Parasitic running. For a fixed mass resolution, 
the resonant search is limited only by the avail¬ 
able statistics. The D*° —>■ D^A' search does not 
require any modifications to the standard LHC run¬ 
ning environment, so it immediately benefits from 
the high data-taking rate (and long run times) al¬ 
ready needed by other LHCb measurements. By 


Strictly speaking, VEPP-3 uses a missing mass strategy. 

One minor difference is that HPS uses strips for its tracking while 
LHCb will use pixels in Run 3. This means that HPS has worse 
hit resolution in the bending plane, so some of the topological 
requirements in App.jB] would not be helpful for HPS. 


contrast, HPS is a dedicated experiment with an 
anticipated runtime of only 3 weeks to cover the 
30-100 MeV mass range. 

• Access to smaller opening angles. As a fixed-target 
experiment, HPS produces A' bosons in the very 
forward direction, effectively within the envelope of 
the beam pipe hole. This means that HPS has no 
sensitivity in the “dead zone” where the A decay 
opening angle is less than 30 mrad. As a colliding 
beam experiment, LHC produces A bosons with 
a sizable transverse momentum kick, such that the 
A trajectory itself goes into the LHCb VELO. This 
allows LHCb to reconstruct much smaller opening 
angles down to around a few mrad, which helps the 
reach at low A masses. 

• Larger Lorentz boosts. The reach in the displaced 
search benehts from large 7 factors (up to a point, 
see below). The median A boost at the LHC 
is roughly three times larger than the maximum 
A boost at HPS. Moreover, there is a tail of A 
events at the LHC which extends to much higher 
boost factors, which can be exploited when com¬ 
bined with the high event rate. 

Note that these last two bullet points are in direct conflict 
at a fixed-target experiment, since going to larger Lorentz 
boosts by using a higher beam energy means that the 
signal has smaller A opening angles. At a colliding beam 
experiment, the opening angle resolution is limited only 
by the hit resolution, so one can in principle exploit larger 
Lorentz boosts up until the point that the displaced A 
signal escapes the detector. From this we conclude that 
HPS is probably close to optimal for a fixed-target dark 
photon search in this mass range. 

IX. SUMMARY 

In this paper, we showed that in Run 3, LHCb can ex¬ 
plore the entire dark photon parameter space between the 
prompt-A' and beam-dump limits for < 100 MeV us¬ 

ing the decay D*° D^A. This reach is possible due 
to the large D*° production rate and sizable A Lorentz 
boost factor at the LHC, combined with the excellent 
vertex/mass resolution and planned triggerless-readout 
system of LHCb. The displaced and resonant strategies 
give complementary coverage of the A parameter space. 
Even if the displaced vertex signature is absent due to 
a modified A lifetime, there is still substantial coverage 
from the resonant search owing to the excellent rUg+g- 
resolution. 

Given the impressive reach below Attid = 142 MeV, 
one might wonder if a similar search could be performed 
at LHCb for dark photons with larger masses. The 
D*o D° A' search relies on the D*° mass constraint to 
suppress backgrounds and to improve the e“'"e“ invariant 
mass resolution. Without these handles, the dark photon 
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search becomes considerably more difficult, but there are 
a number of potential strategies to explore. 

Above the dimuon threshold, a generic search could 
be performed for displaced A' —>• vertices that 

are inconsistent with originating from interactions with 
the detector material. The dimuon invariant mass res¬ 
olution is sufficient without applying a mass constraint, 
and data samples of displaced dimuon vertices are nearly 
background free at LHCb (see, e.g.. Ref. [HI])- A similar 
search was performed by LHCb using Run 1 data that 
looked for the very rare decay —)■ [85]. A key 

challenge for these generic searches, though, is converting 
a cross section bound into a bound on A, since pp col¬ 
lisions provide many potential sources of dark photons 
with uncertain production rates. 

Other radiative charm decays may also offer viable 
A' search channels at LHCb. For example, the yield of 
DO ^ 

7 decays, where both the kaon and pion from 
the K*^—>■ Ktt decay are reconstructed as LONG tracks, 
will be 0(10^°) in Run 3. A search could then be per¬ 
formed for —7- K*^A', with A! —J- e+e”, pApT, and 
7r+7r“. This search is sensitive to mA' ^ 1 GeV and the 
mass constraint would ensure good A' mass resolu¬ 
tion. The A' lifetime is shorter at larger mA'^ however, 
so a displaced search would not be possible over most 
of the mA' range. Nevertheless, such a search is worth 
investigating. Similarly, the decay I?si(2460)“*' —)■ A' 

may provide a viable search for < 500 MeV. 

The common theme for all of the above search strate¬ 
gies is that they benefit greatly from triggerless readout 
with real-time data calibration. These and other searches 
provide strong physics motivation for these upgrades to 
LHCb in Run 3. We encourage more effort in exploring 
how best to exploit these advances in the hunt for dark 
photons and other new particles. 
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Appendix A: Estimating LHCb Tracking 
Performance 


LHCb does not provide public fast simulation software. 
Therefore in this appendix, we estimate the various per¬ 
formance numbers needed for our study from public doc¬ 
uments. 

We start with deriving the resolutions for LONG and 
UP tracks, as defined in Sec. HI A The resolution on Opjp 
for LONG tracks with p > 10 GeV as a function of p 
in the current LHGb detector is provided in Fig. 17 of 
Ref. [Snj. The LONG tracks used in this analysis typi¬ 
cally have momenta in the region where 


« 0.5% 


(Al) 


(n.b. <Tp/p has no strong dependence on p for p < 
10 GeV). The LONG track resolution is expected to im¬ 
prove by about 20% in Run 3, but to be conservative 
we assume the performance will be as measured in the 
current detector. The resolution for UP tracks with mo¬ 
menta typical of electrons coming from A' decays is 


u^P/p « 12%, (A2) 

as can be seen in Fig. 2.5 in Ref. [65] . 

The resolution on the track direction is not explicitly 
published in any LHGb document to our knowledge, but 
we can extract it from published mass resolutions. LHGb 
measures <Jp/p for tracks using various two-body decays, 
e.g. J/A —>■ ■ For the case where p^+ « Pp-, the 

resolution on the track momentum is related to the res¬ 
olution on the J/ijj mass m and decay opening angle a 

by 



\ m J 


\ m J 


(A3) 


which is the same as Eq. (1) of Ref. jSS] but without the 
additional factor of a in the denominator of the rightmost 
term (that is a typo in Ref. [22]). Taking am = 14 MeV 
from Tab. 2, along with ap/p = 0.5% and p « 15 GeV 
from Fig. 17 of Ref. [22], we obtain a a ~ 0.3 mrad, which 
gives a track polar-angle resolution of ag « 0.2 mrad in 
the 9j/^ —>■ 0 limit. 

The value ag = 0.2 mrad should be considered valid 
only for tracks with large momenta. For low-momentum 
tracks, multiple scattering dominates the resolution. We 
simulate a toy model of the Run 3 VELO taking the 
radiation length of the foil and detector modules from 
Ref. [22] ■ The resulting resolution is 

ag K, ^0.2 -I- GeV ^ mrad, a^ « ag cot ff. (A4) 

When the A! is constrained to originate from the 
pp collision, as assumed in deriving Fig. i (7(f) ~ 

( Y ^hit “b ^pp T ) /b ~ 3 mrad (the hit and trans¬ 
verse pp location resolution divided by the radial distance 
to the first hit). 
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We now turn to the resolutions for F-type and P-type 
decays. For two-body F-type decays, both tracks 
are LONG and so one would naively expect dp/p « 
0.5% X ^/2. However, applying the constraint of the 
known mass improves this by about a factor of two 
in our toy simulation. In the —>■ 0 limit, 

~ ~ 0.2mrad, (A5) 

with 

« de cot 0 « 2 mrad. (A6) 

For four-body F-type decays, the resolution is similar 
to the two-body case after applying the mass con¬ 
straint, even if one or two of the tracks are of the UP type. 
For simplicity, the largest dp/p value obtained (when two 
tracks are UP) 

dp ^/p!v 0.5% (A7) 


is used for all F-type decays. 

For P-type decays, the resolution depends strongly 
on the flight distance £. Therefore, improved reso¬ 
lution comes at the expense of signal efficiency. LHCb 
achieves a proper lifetime resolution of dt « 50 fs in 
decays [H5]. Using the relationship 

j + {mef , (A8) 


we can obtain a per-event estimate of the resolution on 
the flight distance. The vertex resolution in the plane 
transverse to the beam line (x — y) is much better than 
along the beam direction (z), e.g. in Ref. [SI] d^jdx ~ 5.5 
for the pp-collision vertex. Therefore, the resolution on 
the transverse flight distance is di.^ ~ cr^sind, and erg « 
(g-fT. ji cos d)(5 .5/-\/2) « tandcr^/d^, where di is obtained 
from Eq. (A8) for each decay. Similarly, cr^ « dijM. 
Both expressions have the expected scaling with flight 
distance. Finally, applying selection criteria that are 
50% efficient gives 


d^^ ~ 2 mrad, d^^ « 30 mrad. (A9) 


In this toy study, we require do < 5 mrad and assume 
additional requirements, e.g. on particle identification, 
are applied that are 90% efficient. In a more sophisticated 
analysis where the full covariance matrix from the vertex 
fit is available, one could choose to cut on dp/p as in 
Ref. [ST]. 

The momentum in P-type decays is solved for as 
described in Sec. lIII Bl For small i it is common that 
is unphysical, i.e. it is larger than the expected to¬ 
tal missing momentum in the rest frame due to the 
resolution on the decay vertex. In such cases, we 
vary the location of the vertex to find the position 
with the smallest vertex that provides a physical 
value. After requiring dg < 5 mrad, we find 


dp^ /p « 9%. 


(AlO) 


It is likely that a more sophisticated selection based on 
dp/p would provide better resolution, though this has 
little impact on the search proposed in this paper. 


Appendix B: Estimating LHCb Vertex Performance 


The resolution on the IP and vertex location are key 
elements of the displaced searches. The one-dimensional 
track IP resolution expected in Run 3 is well approxi¬ 
mated by |M] 


(Tip 


II.O 


I3.IGeV\ 


Pt 




fim. 


(Bl) 


In the 9ji' —>■ 0 limit, this is the IP resolution both in 
and orthogonal to the A' decay plane. We can use dip 
to cross check our estimate of do for tracks. Assuming 
that the location of the pp collision is perfectly measured, 
then d0 « (Tip sin 0/6 mm, where 6 mm is the mean ra¬ 
dial distance of the inner-most track hit [51]. This gives 
de « 2.2GeV/p, which, for typical electron momenta in 
D*o decay s, ag rees to within about 10% with 

our estimate in Eq. (A4) based on mass resolution and a 
toy VELO simulation. 

The resolution on £t in the pre-module search is well 
approximated by 


sin0 PA 9 

cta ~ —— Y^e+ip +^e-ip> 

where o^e+e- is the A' decay opening angle and the A' is 
constrained to originate from the pp collision. The reso¬ 
lution on It: as a function of mA' is shown in Fig. |10[ as 
is the required displacement ndi.^. to achieve « 1 back¬ 
ground event in a given A' signal mass window. The mA' 
dependence is driven by the dependence of (Se+e- on mA' ■ 
In our pre-module displaced search, we require a con¬ 
sistent decay topology. These requirements are as fol¬ 
lows: 

• the A' decay vertex is downstream of the pp colli¬ 
sion; 

• the distance of closest approach between the e+ and 
e~ tracks is consistent with zero; 

• the angle between pA' and the vector formed from 
the pp collision to the A' decay vertex is consistent 
with zero; 

• the IP out of the A' decay plane—defined by the pp 
collision point and the first hits on the and e~ 
tracks—for each electron is consistent with zero. 

In each case, we define consistent with zero as having a 
p-value greater than 1%. Therefore, the efficiency on a 
true displaced A! decay is close to 100%. 

An important point to keep in mind is that virtually all 
fake highly displaced e+e“ vertices that satisfy these con¬ 
sistency requirements are reconstructed with cie+e- much 
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FIG. 10. Distribution of (left) resolution and (right) required displacement nai^ versus mj^i for pre-module A! decays. 


larger than its true value. Therefore, the reconstructed 
value of TOe+e- is larger than its true value in the ab¬ 
sence of substantial bremsstrahlung, and the dominant 
background to the A' signal comes from (more copious) 
e“'"e“ pairs at smaller true invariant mass. 

The opening angle in our search for = 100 MeV is 
on average about the same as for the HPS experiment. 
One can see from Fig. 10 that the pre-module search re¬ 
quirement at 100 MeV is ^ 1 mm or £ > 25 mm. This 
value is similar to the flight distance requirement used 
by HPS to determine the reach of their displaced search, 
leading to a comparable high-side reach at that mass 
value (see Fig.|^. 

For the post-module search, we define IP^ with respect 
to the point where the A! trajectory intersects the 
VETO module, which is the one directly upstream of the 
A' decay vertex location. The resolution on IP^ is given 

by 


CTiP, « + (12^m)2, (B3) 


where is the separation of VETO modules i and 
i + 1, and 12 /rm is the hit resolution. For simplicity, 
we take AZ^ = 25 mm for all modules. In reality, about 
40% of the modules traversed by an A' candidate have 
a larger AZ^. Such cases have slightly better resolution 
since the relative effect of the hit resolution is reduced. 
The resolution on the z position of the A' decay vertex 
in the post-module search is 




1 


Qlg+g— 



+ cr 


2 

e-IPi’ 


(B4) 


where again the A! is constrained to originate from the 
pp collision. The resolution on Zi is shown in Fig. 11 as is 


the required displacement n ■ The consistency require¬ 
ments imposed for the post-module displaced search are 
the same as those used in the pre-module search, with the 
following exceptions: the A! decay vertex is required to 
be downstream of the i**' VETO module, and IP^ is used 
in place of IP for the out-of-decay-plane requirement. 


Appendix C: Belle-II Reach 

The search strategy in this paper can be applied to any 
experiment with a large production rate. As a key 
example, the Belle-B experiment [5H] will collect a data 
sample corresponding to about 50 ab“^ starting in 2017. 
Indeed, Belle-B plans to collect about 10 ab“^ prior to 
the start of Run 3 at the LHC, so they may be able to 
set initial bounds on D*^ —>■ D^A' (or make a discovery) 
prior to LHCb. 

Again neglecting secondary D*° production from b 
hadron decays, the total number of prompt D*^ —>■ 
decays produced at Belle-B will be approximately 

D^-f) « 510pb X 0.381 x 50ab-i«10i°,(Cl) 

where the first number is the inclusive prompt e“''e“ —>■ 
D*^ cross section and the second number is the 
> Z3°7 branching ratio. This is about 500 times less 
than at LHCb in Run 3, though Belle-B will likely be 
able to make use of a larger fraction of decays than 
LHCb. Rescaling the limits derived in this paper, we 
anticipate Belle-B will not be able to probe unexplored 
parameter space using the resonant A’ strategy. 

The main challenge for a Belle-B displaced search is 
that the electrons produced in D*^ —>■ D^A' decays will 
have very low momenta. Using Pythia, we estimate that 
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FIG. 11. Distribution of (left) resolution ctj,. and (right) required displacement ncr^. versus m^/ for post-module A! decays. 


the median value of the lower-momentum electron is only 
60 MeV. This is likely to reduce the tracking efficiency, 
the electron-momentum resolution, and the e"'"e“ vertex 
resolution. The Lorentz 7 factors at Belle-II will also 
be much lower resulting in shorter flight distances. We 


estimate that Belle-II could have sensitivity to displaced 
A decays for small m,A' ■ Given that such sensitivity may 
occur prior to Run 3 at LHCb or elsewhere, we encourage 
a detailed study by the Belle-II collaboration to assess the 
discovery potential. 


[1] L. Okun, Sov. Phys. JETP 56, 502 (1982). 

B166, 473 (1986) 

[2] P. Galison and A. Manohar, Phys. Lett. B136, 279 

[16] S. Gninenko, Phys. Lett. B713, 244 (2012) 

(1984) 

arXiv: 1204.3583 [hep-ph] 

[3] B. Holdom, Phys. Lett. B166, 196 (1986) 

[4] M. Pospelov, A. Ritz, and M. B. Voloshin, Phys. Lett. 

[17] F. Archilli et al. (KLOE-2 collaboration), Phys. Lett. 
B706, 251 (2012), arXiv:1110.0411 [hep-ex] 

B662, 53 (2008), arXiv:0711.4866 [hep-ph] 

[18] D. Babusci et al. (KLOE-2 collaboration), Phys. Lett. 

[5] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, 

B720, 111 (2013), arXiv: 1210.3927 [hep-ex] 

and N. Weiner, Phys. Rev. D79, 015014 (2009) 

[19] R. Essig et al., in Community Summer Study 2013: 

arXiv:0810.0713 [hep-ph] 

Snowmass on the Mississippi (CSS2013) Minneapo¬ 

[6] J. D. Bjorken, R. Essig, P. Schuster, and N. Toro, Phys. 

lis, MN, USA, July 29-August 6, 2013 (2013) 

Rev. D80, 075018 (2009), arXiv:0906.0580 [hep-ph] 

arXiv:1311.0029 [hep-ph] 

[7[ G. Bernardi, G. Carugno, J. Chauveau, F. Dicarlo, 

[20] B. Aubert et al. (BaBar collaboration), in Proceedings, 

M. Dris, et al, Phys. Lett. B166, 479 (1986) 

34th International Conference on High Energy Physics 

[8] R. Meijer Drees et al. (SINDRUM I collaboration), Phys. 

(ICHEP 2008) (2008) arXiv:0808.0017 [hep-ex] 

Rev. Lett. 68, 3845 (1992) 

[21] B. Batell, M. Pospelov, and A. Ritz, Phys. Rev. D80, 

[9[ P. Astier et al. (NOMAD collaboration), Phys. Lett. 

095024 (2009), arXiv:0906.5614 [hep-ph] 

B506, 27 (2001), arXiv:hep-ex/0101041 [hep-ex] 

[22] P. deNiverville, M. Pospelov, and A. Ritz, Phys. Rev. 

[10] S. Gninenko, Phys. Rev. D85, 055027 (2012) 

D84, 075020 (2011), arXiv:1107.4580 [hep-ph] 

arXiv:1112.5438 [hep-ph] 

[11] P. Adlarson et al. (WASA-at-COSY collaboration), Phys. 

[23] B. Wojtsekhowski, D. Nikolenko, and 1. Rachek, (2012), 
arXiv: 1207.5089 [hep-ex] 

Lett. B726, 187 (2013), arXiv: 1304.0671 [hep-ex] 

[24] Y. Kahn and J. Thaler, Phys. Rev. D86, 115012 (2012) 

[12] G. Agakishiev et al. (HADES collaboration), Phys. Lett. 
B731, 265 (2014), arXiv:1311.0216 [hep-ex] 

arXiv: 1209.0777 [hep-ph] 

[25] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, Phys. 

[13[ A. Adare et al. (PHENIX collaboration), Phys. Rev. 
C91, 031901 (2015), arXiv: 1409.0851 [nucl-ex] 

Rev. D88, 114015 (2013), arXiv: 1307.6554 [hep-ph] 

[26] B. Batell, R. Essig, and Z. Surujon, Phys. Rev. Lett. 

[14[ J. Batley et al. (NA48/2 collaboration), (2015), 

arXiv:1504.00607 [hep-ex] 

113, 171802 (2014), arXiv: 1406.2698 [hep-ph] 

[27] Y. Kahn, G. Krnjaic, J. Thaler, and M. Toups, Phys. 

[15] F. Bergsma et al. (CHARM collaboration), Phys. Lett. 

Rev. D91, 055006 (2015), arXiv: 1411.1055 [hep-ph] 











19 


[28] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, Phys. 
Rev. D91, 094026 (2015), arXiv:1411.1404 [hep-ph] 

[29] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, 
(2015), arXiv:1505.00011 [hep-ph] 

[30] H. Merkel et al. (A1 collaboration), Phys. Rev. Lett. 106, 
251802 (2011), arXiv:1101.4091 [nucl-ex] 

[31] H. Merkel et al, Phys. Rev. Lett. 112, 221802 (2014) 
arXiv:1404.5502 [hep-ex] 

[32] S. Abrahamyan et al. (APEX collaboration), Phys. Rev. 
Lett. 107, 191804 (2011), arXiv: 1108.2750 [hep-ex] 

[33] A. Konaka, K. Imai, H. Kobayashi, A. Masaike, 
K. Miyake, et al., Phys. Rev. Lett. 57, 659 (1986) 

[34] E. Riordan, M. Krasny, K. Lang, P. De Barbaro, 
A. Bodek, et al., Phys. Rev. Lett. 59, 755 (1987) 

[35] J. Bjorken, S. Ecklund, W. Nelson, A. Abashian, 
C. Church, et al., Phys. Rev. D38, 3375 (1988) 

[36] A. Bross, M. Crisler, S. H. Pordes, J. Volk, S. Errede, 
et al., Phys. Rev. Lett. 67, 2942 (1991) 

[37] M. Davier and H. Nguyen Ngoc, Phys. Lett. B229, 150 
(1989) 

[38] C. Athanassopoulos et al. (LSND collaboration), Phys. 
Rev. C58, 2489 (1998), arXiv:nucl-ex/9706006 [nucl-ex] 

[39] S. Adler et al. (E787 collaboration), Phys. Rev. D70, 
037102 (2004), arXiv:hep-ex/0403034 [hep-ex] 

[40] A. V. Artamonov et al. (BNL-E949 collaboration), Phys. 
Rev. D79, 092004 (2009), arXiv:0903.0030 [hep-ex] 

[41] R. Essig, R. Harnik, J. Kaplan, and N. Toro, Phys. Rev. 
D82, 113008 (2010), arXiv: 1008.0636 [hep-ph] 

[42] J. Blumlein and J. Brunner, Phys. Lett. B701, 155 
(2011), arXiv:1104.2747 [hep-ex] 

[43] J. Blumlein and J. Brunner, Phys. Lett. B731, 320 
(2014), arXiv:1311.3870 [hep-ph] 

[44] D. Hanneke, S. Fogwell, and G. Gabrielse, Phys. Rev. 
Lett. 100, 120801 (2008), arXiv:0801.1134 [physics.atom- 
ph] 

[45] G. F. Giudice, P. Paradisi, and M. Passera, JHEP 11, 
113 (2012), arXiv: 1208.6583 [hep-ph] 

[46] M. Endo, K. Hamaguchi, and G. Mishima, Phys. Rev. 
D86, 095029 (2012), arXiv: 1209.2558 [hep-ph] 

[47] M. Pospelov, Phys. Rev. D80, 095002 (2009) 

arXiv:0811.1030 [hep-ph] 

[48] J. B. Dent, F. Ferrer, and L. M. Krauss, (2012), 
arXiv:1201.2683 [astro-ph.CO] 

[49] D. Kazanas, R. N. Mohapatra, S. Nussinov, V. L. 
Teplitz, and Y. Zhang, Nucl. Phys. B890, 17 (2014) 
arXiv:1410.0221 [hep-ph] 

[50] K. A. Olive et al. (Particle Data Group), Chin. Phys. 
C38, 090001 (2014) 

[51] “LHCb Trigger and Online Technical Design Report,” 
(2014), LHCb-TDR-016. 

[52] O. Moreno, (2013), arXiv:1310.2060 [physics.ins-det] 

[53] T. Sjostrand et al., Comput. Phys. Commun. 191, 159 
(2015), arXiv:1410.3012 [hep-ph[ 

[54] R. Aaij et al. (LHCb collaboration), Nucl. Phys. B871, 
1 (2013), arXiv:1302.2864 [hep-ex] 

[55] R. Aaij et al. (LHCb collaboration), “LHCB-PAPER- 
2015-041, in preparation,”. 

[56] G. A. Miller and P. Singer, Phys. Rev. D37, 2564 (1988) 

[57] 1. W. Stewart, Nucl. Phys. B529, 62 (1998), arXiv:hep- 
ph/9803227 [hep-ph] 


[58] A. A. Alves Jr. et al. (LHCb collaboration), JINST 3, 
S08005 (2008) 

[59] R. Aaij et al. (LHCb collaboration), Int. J. Mod. Phys. 
A30, 1530022 (2015), arXiv: 1412.6352 [hep-ex] 

[60] R. Aaij et al, JINST 9, P09007 (2014), arXiv:1405.7808 
[physics.ins-det] 

[61] M. Adinolfi et al., Eur. Phys. J. C73, 2431 (2013) 
arXiv:1211.6759 [physics.ins-det] 

[62] A. A. Alves Jr. et al, JINST 8, P02022 (2013) 
arXiv:1211.1346 [physics.ins-det] 

[63] R. Arink et al., JINST 9, P01002 (2014), arXiv:1311.3893 
[physics.ins-det] 

[64] “LHCb VELO Upgrade Technical Design Report,” 
(2013), LHCb-TDR-013. 

[65] “LHCb Tracker Upgrade Technical Design Report,” 
(2014), LHCb-TDR-015. 

[66] R. Aaij et al. (LHCb collaboration), JHEP 07, 014 
(2014), arXiv:1405.2797 [hep-ex] 

[67] R. a. Aaij (LHCb collaboration), JHEP 10, 005 (2014) 
arXiv: 1408.1299 [hep-ex] 

[68] H. W. Koch and J. W. Motz, Rev. Mod. Phys. 31, 920 
(1959) 

[69] “LHCb PID Upgrade Technical Design Report,” (2013), 
LHCb-TDR-014. 

[70] HPS collaboration, “Proposal for 2014-2015 run,”. 

[71] A. Borsellino, Phys. Rev. 89, 1023 (1953) 

[72] W. D. Hulsbergen, Nucl. Instrum. Meth. A552, 566 
(2005), arXiv:physics/0503191 [physics] 

[73] R. Aaij et al. (LHCb collaboration), JHEP 04, 064 
(2015), arXiv:1501.03038 [hep-ex] 

[74] S. Stone, “private communication,”. 

[75] LAPPD collaboration, “http://psec.uchicago.edu,”. 

[76] M. Williams, JINST 10, P06002 (2015) 

arXiv: 1503.04767 [hep-ex] 

[77] R. Essig, P. Schuster, N. Toro, and B. Wojtsekhowski, 
JHEP 02, 009 (2011), arXiv: 1001.2557 [hep-ph] 

[78] T. Beranek, H. Merkel, and M. Vanderhaeghen, Phys. 
Rev. D88, 015032 (2013), arXiv: 1303.2540 [hep-ph] 

[79] M. Freytsis, G. Ovanesyan, and J. Thaler, JHEP 01, 
111 (2010), arXiv:0909.2862 [hep-ph] 

[80] J. Balewski, J. Bernauer, W. Bertozzi, J. Bessuille, 
B. Buck, et al., (2013), arXiv: 1307.4432 

[81] B. Echenard, R. Essig, and Y.-M. Zhong, JHEP 01, 113 
(2015), arXiv:1411.1770 [hep-ph] 

[82] S. Alekhin et al., (2015), arXiv: 1504.04855 [hep-ph] 

[83] S. Gardner, R. J. Holt, and A. S. Tadepalli, (2015), 
arXiv: 1509.00050 [hep-ph] 

[84] R. Aaij et al. (LHGb collaboration), (2015), 
arXiv: 1508.04094 [hep-ex] 

[85] R. Aaij et al. (LHCb collaboration), JHEP 01, 090 
(2012), arXiv:1209.4029 [hep-ex] 

[86] R. Aaij et al. (LHCb collaboration), Phys. Rev. Lett. 
110, 101802 (2013), arXiv:1211.1230 [hep-ex] 

[87] R. Aaij et al. (LHCb collaboration). Nature Phys. 11, 17 
(2015), arXiv:1504.01568 [hep-ex] 

[88] Belle collaboration, “Belle 11 technical design report,” 
KEK Report 2010-1. 





